Charge and Magnetization Inhomogeneities in Diluted Magnetic Semiconductors 
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It is predicted that III-V diluted magnetic semiconductors can exhibit stripe-like modulations of magnetization 
and carrier concentration. This inhomogeneity results from the strong dependence of the magnetization on the 
carrier concentration. Within Landau theory, a characteristic temperature T* below the Curie temperature is 
found so that below T* the equilibrium magnetization shows modulations, which are strongly anharmonic. 
Wavelength and amplitude of the modulation rise for decreasing temperature, starting from zero at T* . Above 
T* the equilibrium state is homogeneous, but the coupling between charge and magnetization leads to the 
appearance of an electrically charged layer in domain walls. 

PACS numbers: 75.50.Pp, 75.30.Fv, 75.60.Ch, 75.10.Hk 
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Introduction. — Diluted magnetic semiconductors (DMS) 
are investigated extensively as promising materials for spin- 
tronics applications 1 1 , 2] and because of their unique physical 
properties |3j,|4j]. Since the magnetic interaction is mediated 
by the carriers, the magnetization and the Curie temperature 
increase for increasing carrier concentration j^.lril7il8ll9[ floll . 
In fact, the magnetization can be changed in situ by tuning the 
carrier concentration with a gate voltage 

BE!. 

Coupling between magnetism and carrier concentration can 
lead to inhomogeneous equilibrium states, as found in man- 
gam tes , nickelates and cuprates in. The present 
Letter analyzes this possibility for III-V DMS, concentrat- 
ing on stripe-like, one-dimensional variations of the magne- 
tizations. We employ a Landau theory for the coupled mag- 
netic and charge degrees of freedom. This approach is valid 
on length scales on which the impurity distribution can be 
treated as homogeneous The characteristic length scale 
is n^J 3 , where riyi n is the density of Mn impurities fil lT^I . 
Finally, we discuss possible experiments. 

Magnetic domains have also been observed in (Ga,Mn)As 
flcl Mm . They are formed to reduce the dipolar energy, 
as in other ferromagnets. One can expect the coupling be- 
tween magnetization and carriers to lead to an inhomogeneous 
charge profile in a domain wall. This question is addressed in 
the final part of this Letter. 

Landau theory. — We write the Hamiltonian as a functional 
of magnetization m and deviation of carrier density from its 
spatial average, Sn = n — n. Charge neutrality requires 
/ d 3 r Sn = 0. The magnetic part has the usual form 
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where di = d/dri and summation over i is implied. The 
mean-field Curie temperature is determined by a = 0. Since 
experimentally the Curie temperature depends approximately 
linearly on carrier concentration, we expand a = a' (T — T c — 
n Sn), where T c is the Curie temperature for Sn = 0. This 
dependence of a provides the coupling between magnetism 
and carrier concentration in our model and is responsible for 
the physics discussed in the following. Since the equilibrium 



magnetization for constant a is mo = -J —a/ [3, larger magne- 
tization is favored in regions with higher carrier concentration 
(note a < in the ferromagnetic phase). 

The second ingredient for our model is the screened 
Coulomb energy due to the charge inhomogeneity, 
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The total Hamiltonian is H = H m + H$ n . 

We first discuss qualitatively what kind of equilibrium 
states we expect from H. Any inhomogeneous charge dis- 
tribution increases Hs n . On the other hand, the contribution 
from H rn is not obvious, since the first term is negative for 
T < T c + n Sn. We will see that the magnetic energy de- 
crease in regions of higher carrier concentration and magne- 
tization can outweigh the increase in regions of lower Sn and 
to and even the increase in electrostatic energy. In that case, 
the equilibrium state is indeed inhomogeneous. We consider 
stripe-like, one-dimensional modulations. Two- and three- 
dimensional patterns seem less likely, because they contain 
more regions with large magnetization gradients for a given 
inhomogeneity length scale, which increase the energy due to 
the gradient term in H m . One could expect the inhomogeneity 
to take the form of stripe domains 1 16] with alternating mag- 
netization. However, we will see that the equilibrium solution 
shows a magnetization modulation without sign change. 

We now turn to the formal derivation of equilibrium states. 
It is convenient to express Hs n in terms of the electrostatic 



potential 4>. With (A - 
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DMS) we obtain the total Hamiltonian 
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Equilibrium configurations are given by minima of H subject 
to the constraint of charge neutrality, J d 3 r (A — rQ 2 )cf> = 0, 
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FIG. 1: (color online). Phase diagram for periodic solutions for the 
magnetization in terms of m m i n and m* = m(0) for m S m g = 1. 



Type 1 solutions exist for m* 
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and be- 



come a special case of type 2 for m* < m s j ng . Type 2 solu- 
tions exist for m* < min(m m i n , m s i ng ) and type 3 solutions for 
7n s i n g < m* < (2m^ nin — f7igi n g) ly/2 . The homogeneous solutions 
are also shown. Inset: Schematic plot of the denominator in Eq. jTJ 
showing the values assumed by m for type 2 and 3 solutions. 



which is implemented with a Lagrange multiplier. Introducing 
the averaged squared magnetization to 2 and the rescaled po- 
tential $ = (a'r]€oe/e) <j> we obtain two coupled Euler equa- 
tions for m and $. Eliminating <f> from the first we find 
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with a = a' 2 i] 2 eoe/e 2 , b = 7, c = a' (T — T c ), d = (3. Equa- 
tion (0J retains information about the coupling to the carrier 
density since a oc rj 2 . The only bounded solution of Eq. (|5} 
is $ = (a/2) (to 2 — to 2 ). The equations support the homo- 
geneous mean-field solution to 2 = to 2 , = — c/d > 0, $ = 
for T < T c . While Eqs. J4l5i contain five parameters, it is 
sufficient to vary only two to obtain all possible solutions up 
to rescaling. We choose a/dr 2 , oc rj 2 and cr^/b oc T — T c , 

Periodic solutions. — For periodic, collinear solutions that 
only depend on x, Eq. Q becomes 
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This is an integro-differential equation due to the term to 2 = 
(1/A) j^dxm 2 (x) in c', where A is the wavelength, c' 
is determined selfconsistently below. Boundary conditions 
m(0) = to* and d x m(Q) = are imposed, where to* will 
be obtained by minimizing the energy. 

Using standard methods, we obtain the explicit integral for 
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FIG. 2: (color online), (a) Phase diagram for periodic solutions for 
the magnetization in terms of cr^/b oc T — T c and initial magnetiza- 
tion to* = to(0) for ajdr\ = 1. The symbols are as in Fig.Q In the 
distorted triangle with corners A,B,0, two solutions with different 
wavelength coexist, (b) Phase diagram for periodic magnetization 
solutions in terms of crl/b (x T -T c and drl/ a oc m„ in „. The 
various solutions exist to the left of the respectively lines. 



the inverse function on the interval [—A/2, A/2], 
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This expression satisfies d x m(0) — 0, since dx/dm diverges 
at to = to*. To obtain a periodic function, the denominator 
must have another zero at the next extremum of m(x) at x — 
±A/2. Beyond ±A/2 the solution continues periodically. 

A special role is played by the magnetization value m s i ng = 
\Jbja: Here dx/dm vanishes so that coordinates x beyond 
this point normally cannot be reached and there is no solu- 
tion for all x. However, a solution (here called type 1) cross- 
ing to = TOsing is possible if numerator and denominator 
vanish simultaneously. This solution is m(x) = (2m 2 nin — 
to 2 ^) 1 / 2 cos(-y/ d' /2a x) where m m j n = yj—c'/d' is the 
minimum of the denominator. 

For all other periodic solutions to must be either larger 
or smaller than TO s ; ng everywhere. For to* < TO s i ng , pe- 
riodic solutions (type 2) oscillating between to* and —to* 
with zero average exist if to* < TO m ; n . For to* > m s i ng , 
periodic solutions {type 3) oscillating around m m \ n exist for 
c'm* 2 + d'm* 4 /2 < c'm 2 ing + d'mf ing /2, see Fig.[T] 

Next, the phase diagram is mapped back onto the param- 
eters of the Euler equation ©, determining to 2 selfconsis- 
tently. The results are shown in Fig.|2ja) for a/dr^ = 1. The 
diagram for other values has the same topology. Figure |2jb) 
shows the resulting phase diagram in terms of cr\/b and gen- 
eral dr 2 /a oc TO 2 ing . The lines A, B, C show the shift of the 
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crossing points marked A, B, C in Fig.[2ja) with dr^/a. Type 
3 solutions exist to the left of point A, i.e., for mo > m s i n g- 
From Eqs. (|3}-l|3} we find the average energy density 
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For the homogeneous solution we obtain the standard result 
e = £hom = —c 2 /4d. Numerical evaluation shows that type 
1 and 2 (type 3) solutions always have higher (lower) energy 
than the homogeneous solution. Among type 3 solutions the 
energy is minimized by the maximum amplitude, where m 
comes arbitrarily close to m s i n g- 

The mean-field magnetization of our DMS model is thus 
zero for T > T c , homogeneous for T* <T < T c , where 



T* = T c 
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corresponds to line A in Fig.|2jb), and a periodic spin-density 
and charge-density wave for T < T* . The dipolar interaction 
omitted here favors m lying in the yz plane. The magnetiza- 
tion and potential show sharp cusps at the minima of m, see 
the inset in Fig. [3] The cusps lead to negative peaks in the car- 
rier density, which become 5-functions for m* — > m s i ng . This 
divergence is cut off by the condition of non-negative hole 
concentration. Since the amplitude, wavelength, and energy 
approach finite values for m* — > TO S j ng , the Landau theory 
gives a good impression of the profile, except for some broad- 
ening of the cusps. 

The optimum solution can be written down explicitly, 



m(x) = y2ni 2 . 
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for < x < A and periodically repeated. Here, (2m 
m 2 ;^) 1 / 2 sin 8 = m S j ng . From m(x) one can obtain expres- 
sions for the wavelength A = 2(2a/d') 1/2 (tt/2 - 6), the av- 
erage magnetization m, the peak-to-peak amplitude <5m pp — 
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and the energy, see Fig. [5] Note 



that m is nonzero for all T < T c . Close to T* the wavelength 
becomes small. In this regime, the continuum theory breaks 
down, since A is not large compared to the disorder length 
scale. Figure [3] also shows that the fundamental length scale 
is the screening length tq- The inset shows a typical solution. 

It is important to check whether the periodic solution can 
occur in real DMS. For that, T c — T* should be small. Equa- 
tion ij^Ji shows that this is the case for high dielectric constant, 
small spin stiffness, strong dependence of T c on carrier con- 
centration, and rapid onset of magnetization below T c . For 
(Ga,Mn)As we estimate T* by comparing experimentsO, 0] 
to mean-field theory for homogeneous magnetization 1 14|] and 
to spin-wave theory GJ]. We find T c - T* of the order of 
10 K. The properties of (Ga,Mn)As vary strongly with Mn 
concentration and growth procedures. In particular, T c — T* 
is inversely proportional to the square of the shift of T c with 
carrier concentration, ?y 2 . In Ref. |9], r\ ~ 5.4 x 10 5 KA 3 , 
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FIG. 3: (color online). Wave length A, average magnetization m, 
peak-to-peak amplitude Sm pp , and gain in energy density e — e^om 
for the periodic magnetization with lowest energy. The magnetiza- 
tion for the homogeneous solution is also shown. The unit of mag- 
netization is m s = yjbfd/ro. a/dr^ = 1 is assumed. Inset: mag- 
netization and excess carrier density for cr^/b = —4. 



which was used for the above estimate, whereas in Ref. 0] 
t] ~ 1.5 X 10 5 KA 3 , which would increase T c - T*. 

Figure[3]suggests that measurements of the average magne- 
tization, which have been performed extensively, are unlikely 
to find evidence for the inhomogeneous state. For that, probes 
sensitive to the spatial variation are required. For example, 
the magnetic modulation should be observable in neutron- 
scattering experiments. In real space, magnetic scanning- 
tunneling microscopy (STM) and, for large A, scanning Hall 
probe experiments flrill or magneto-optical techniques ifTvtl 
are promising. Conversely, the modulation in carrier con- 
centration should be observable in optical reflection or trans- 
mission for large enough A. It also leads to a modulation of 
the local density of states which could be probed by STM. 
The smoking gun experiment would be to look for charge and 
magnetization modulations of the same wavelength. 

Domain walls. — Finally, we study the effect of magnetiza- 
tion-carrier coupling on domain walls HI <tI 1 1 '711 - We restrict 
ourselves to solutions that are homogeneous in the y, z direc- 
tions. Equations ( 1415 \ are solved under the boundary condi- 
tions lirxi x _,±oo m(x) = ±moz, where z is the unit vector in 
the z direction. Since m 2 (x) only deviates appreciably from 



m,Q in a finite interval, we have m 2 



-c/d in the limit 



of infinite system size, L — > oo. However, it turns out that 
charge neutrality can only be satisfied by keeping terms of or- 
der 1/L in m 2 . One such term comes from the region far from 
the wall, where we write m(L/2) = mo + mi/L. This means 
that the enhanced carrier density compensating the reduction 
in the wall is spread out over the bulk. 

With the additional condition m(0) = we obtain 



dm 



d'm A /2 



cmf. 



d'ml/2 



(11) 



4 




FIG. 4: (color online). Magnetization and electrostatic potential for 
typical domain- wall solutions for ajdr\ — 1. The curves are for 
crl/b = -0.1, -0.3, -0.5, -0.7, -0.9, -0.99. The asymptotical 
values mo = m(x — > oo) are indicated by triangles. Inset: Energy 
density Act of a domain wall as a function of coupling strength for 
the same values of c. 



tial decrease is dominated by the 1/L term in Am 2 (x) far 
from the wall, i.e., from the redistribution of carriers. In this 
regime, domain walls are (slightly) less costly than they would 
be without coupling. The strong increase mostly comes from 
the increased width due to Coulomb repulsion. Domain walls 
could be observed in the real-space experiments discussed 
above. The charged layer should also affect electronic trans- 
port through domain walls. 

Conclusions. — The carrier-concentration dependence of 
the magnetization in DMS introduces a characteristic tem- 
perature T* < T c such that the mean-field magnetization rn 
and excess carrier density 5n show periodic modulations for 
T < T*, whereas m is homogeneous and 5n = above T*. 
T c - T* can be of the order of 10 K in p-type DMS. The mod- 
ulation is strongly anharmonic, and amplitude and wavelength 
increase for decreasing temperature, starting from zero at T*. 
For T > T* the equilibrium state is homogeneous, but the 
coupling between magnetism and carrier concentration leads 
to the appearance of a negatively charged layer in the vicinity 
of a domain wall for p-type DMS. 



withe' = c—{a/2rf j )(c/d) and d! = d—a/2r 2 . The integrand 
must be free of singularities for < m < mo, which implies 
?7io < msjjw. Thus domain-wall solutions only exist for T* < 
T <T C 03, see Fig.EJa). This indicates that T c - T* is large 
for samples which show domains at low temperatures. 

Equation can be integrated explicitly. It also yields 
an expression for the typical width t; w = nio/d x m(0) of a 
domain wall, = -2bd/c(d - a/2r 2 ). £ w increases for 
increasing coupling a oc rj 2 between magnetism and carriers 
due to their Coulomb repulsion. Figure |4] shows m{x) and 
for typical domain-wall solutions. The excess carrier 
concentration Sn oc (r^ 2 — A)$ is negative in the domain 
wall, where the magnetization is reduced. 

The areal energy density of the domain wall is obtained by 
integrating the energy density over x, where corrections to 
m(x) of order 1/L are again relevant, 
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with Am 2 = m 2 (±oo) — m 2 (x). The dependence of Act 
on the coupling a oc rj 2 is shown in the inset in Fig. |4] Act 
first decreases with increasing coupling and then increases 
again, finally diverging as d! = d — a/2r^ goes to zero. For 
larger | cr§ / | the divergence is not reached, since the condi- 
tion mo = m S j ng (i.e., T — T*) is satisfied first. The ini- 
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